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ABSTRACT 

VLA surface brightness and spectral index maps of the evolving extended emission of 
the triple symbiotic star CH Cygni are presented. These are derived from observations 
at 4.8, 8.4 and 14 GHz between 1985 and 1999. The maps are dominated by thermal 
emission around the centrai bright peak of the nebula, but we also find unambiguous 
non-thermal emission associated with the extended regions. Our observations confimi 
that this is a jet. The centrai region has been associated with the stellar componcnts 
through HST imaging. If the jets are due to ejection events at outburst, expansion 
velocities are consistent with those from other measurement methods. We propose that 
the non-thermal emission is caused by material ejected in the bipolar jets interacting 
with the circumstellar wind envelope. The resulting shocks lead to locai enhancements 
in the magnetic field from the compact componcnt of the order 3 mG. 

Key words: binaries: symbiotic - circumstellar matter - stars: imaging - stars: 
individuai: CH Cygni - stars: late-type - radio continuum : stars 



1 INTRODUCTION 

Symbiotic stars occupy an extreme and relatively poorly 
understood region of the binary classifìcation scheme. The 
name was coined by Merril (1941) to describe stars which 
appeared to have a combination spectra: that of high ex- 
citation lines usually associated with a hot, ionized nebula 
superimposed on a cool continuum with prominent absorp- 
tion features consistent with a late-type star. At present 
they are understood to be interacting binaries (with orbitai 
periods of a few to tens of years) consisting of a cool giant 
losing material mostly via the stellar wind and a hot, lu- 
minous compact object which ionises a portion of the cool 
component wind (Seaquist et al., 1984, hereafter the STB 
model). Such a state of affairs represents the so-called qui- 
escent phase, which can be interrupted by periods of activity. 
The active phases start with an eruption of the hot star, an 
event indicated photometrically by an increase of the star's 
brightness by 2-6 mag, and/or spectroscopically by high ve- 
locity broad emission features from the centrai star. Both 
radio and HST imaging can directly resolve the remnants 
of such dramatic events (Hack & Parsece, 1993; Eyres et al. 
1995; Kenny et al. 1996; Richards et al., 1999; Eyres et al., 
2001; Eyres et al. in preparation, hereafter Paper II). 

The symbiotic star CH Cygni displays particularly 



complex behaviour. Optical spectroscopic studies by 
Mikolajewski et al (1987) showed that the orbits of the stars 
within the system are likely to be coplanar and eclipsing, 
with eclipses separated by around 5700 days. Later studies 
(Mikolajewski et al, 1990 and references therein) confirmed 
this period. 

Further spectral studies of the system (Hinkle et al. 
1993) suggested that CH Cyg is probably a triple-star sys- 
tem consisting of an inner 756-day period binary which is 
orbited by an unseen G-K dwarf on a 5300 day orbit. Further 
photometric studies (Skopal et al., 1996) discovered short, 
756 day, period eclipses which show that ali three stars in 
the system are likely to be in coplanar orbits. 

The activity of CH Cyg is seen best in the U light- 
curve (Skopal, 1998), which indicates a brightening from 
U ~11-12 to 5 (Fig. [l]). The drop in the optical magni- 
tude of CH Cygni at the times of the long-period eclipses 
have been marked on this figure. Each outburst was accom- 
panied by high velocity broad emission features consistent 
with mass outflows. During 1984-85 the material was ejected 
at ~600-2500kms _1 (Mikolajewski & Tomov 1986), whilst 
the 1992-95 active phase was characterized by sporadic and 
in part bipolar outflow at ~ 1000- 1600 km s -1 (Leedjàrv & 
Mikolajewski, 1996; Skopal et al. 1996; Ijima 1996) and, fi- 
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Figure 1. The U light-curve for CH Cyg covering the period 1967 — 2000. The eclipses in the long, 14.5-year, period orbit are marked 
by E. The timings of our VLA observations pertinent to this paper are denoted by dashed bars. 



Table 1. VLA Observations of CH Cygni. 



Date 


Band 1 


Band 2 


Time on Source 


Map lnoisc^j 


Time on Source 


Map 2 noise 




(GHz) 


(GHz) 


at Band l(s) 


(mjy beam -1 ) 


at Band 2(s) 


(mjy beam -1 ) 


1986 Mar 20 


4.8 


15 


2290 


0.053 


9930 


0.083 


1988 Oct 16 


4.8 


15 


2600 


0.032 


5220 


0.073 


1995 Aug 20 


1.8 


8.5 


2550 


0.021 


5620 


0.063 


1999 Sep 26 


4.8 


8.5 


3210 


0.077 


2400 


0.133 



referred to as (T 



nally, during the recent, 1998-2000, outburst mass outflows 
at about 1000 km s -1 were observed (Paper II). The out- 
flows may be correlated with a significant increase of the ra- 
dio emission and the radio light curves during these periods 
correlate well with the optical ones (Kenny et al. 1996). The 
1984 mass ejection has been linked to emergence of bipolar 
emission which could be from either an expanding torus or 
collimated jets (Taylor et al. 1986) . Later radio observations 
revealed a complex morphology, although stili exhibiting a 
bipolar structure (Kenny et al. 1996). 

In this paper we present radio surface brightness and 
spectral index maps taken from 15 years of observations of 
CH Cygni by the Very Large Array. 



2 OBSERVATIONS 

CH Cygni has been observed in the highest resolution "A" 
configuration several times with the VLA since 1984 and 
multiple frequency observations were made on many of these 
occasions, four of which are detailed in Table |l|. 

Data reduction was carried out according to standard 
procedures using the AIPS software package (Greisen, 1999). 
The fiux density of the phase reference source was deter- 
mined, and the data calibrated to correct phase and am- 



plitude instabilities. These corrections were applied to the 
observations of CH Cygni and a map produced by deconvolv- 
ing the instrumentai response (the so-called "dirty beam") 
from the map using the well-established CLEAN algorithm 
(Hogbom, 1974). Due to the low radio fiux density of the 
source no self-calibration was possible. The maps were con- 
volved with a beam size of 0.5x0.5 arcsec at ali observing 
bands. This size was the mean of the the major and mi- 
nor axes of the beam at 5GHz, in order that comparison of 
the size and shape of features at different epochs would be 
meaningful. 



Observations of the same object at different frequencies 
have a different u, v piane coverage. This leads to features 
of different spatial scales being picked out at the different 
frequencies, causing problems when taking a ratio between 
two maps. In order to avoid this, only the u, v range cov- 
ered by both sets of data was used to create the maps. This 
has the undesirable effect of decreasing the signal-to-noise 
ratio but does mean that the maps can be compared. The 
AIPS task "comb" was then used to produce a third map 
(at a resolution of 0.17 arcseconds per pixel) which traces 
the pixel-by-pixel spectral index variation. We use the con- 
vention that the radio fiux at frequency v, S v oc v a so that 
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Figure 2. HST WFPC2 image of CH Cygni from 1999 August 13. 
The point spread function of the centrai star has been subtracted 
from this image for clarity. 



a positive spectral index (a) represents an increase in flux 
with frequency. 

Hubble Space Telescope imaging was used to place the 
stellar components on the radio map. The observations were 
made as part of a GO programme 8330 on symbiotic stars. 
Three orbits were allocated to CH Cygni. A WFPC2 image 
(Fig. §) taken through the F502N filter (Biretta 1996) on 
1999 August 13 was aligned to the VLA map of 1999 Septem- 
ber. The best agreement was achieved when the radio peak 
coincided with the star in the optical domain (this is dis- 
cussed in more detail in Paper II). 



3 RESULTS 

The radio brightness and spectral index maps are shown 
in Fig. |^. The brightness maps are dominated by a centrai 
bright peak surrounded by an extended region. The mor- 
phology of these extensions varies from one epoch to the 
next. In particular, after the initial bipolar extension ob- 
served in 1986 the extension appears to be more pronounced 
on the SE side of the centrai peak. In addition to this the 
position angle of the extension has changed (Crocker et al, 
in preparation) . 

In each spectral index map, there is a wide range of in- 
dices visible. In particular, the 1986 map (Fig. ^, top) shows 
a distinct trend with positive a in the centrai region and 
negative a in the outlying areas. This strongly suggests the 
existence of thermal emission around the stellar components 
with extended non-thermal emission from the well-resolved 
knots. 

Cross-sections were taken through the peak of each 
spectral index map, along the axis of maximum elongation 
in the 5GHz map. In order to reduce the uncertainties, a sec- 
tion was taken on each side of the initial one, giving three 
parallel sets of data. At each point along the slice, the three 
pixels were combined to find a weighted mean and then used 
to create the final cross-section. The weighted mean cross- 
sections are shown in Fig. ^, where the positions are given 
relative to the centrai bright peak in the 8.5GHz or 15GHz 
radio maps, whichever was available. A feature common to 



ali the epochs is a centrai peak, accompanied by wings of 
lower spectral index. 

The first plot, from 1986, has a peak value of a — 0.42± 
0.03 which is consistent with that expected from a Constant 
outfiow ionized wind, assuming a variation in density of p oc 



-0 



with j3 



The decrease in a with distance from 



the centre favours a geometry which is optically thick along 
lines of sight that pass through a dense core. In the extended 
regions, the mean value of a is —0.23 ± 0.03. 

By 1988, the spectral index of the centrai region had 
increased to a = 1.02 ± 0.04 as if /3 > 2, which can be 
attributed to the photoionised nebula dominating over the 
wind. Such a value is predicted by the STB model for 
lower ionisation parameters. No extended emission is seen 
at 15GHz. This observation was made at a time when the 
optical activity was very low (Fig. |l|) and no outburst had 
been observed for some time. 

The 1995 cross-section has a peak of similar width to the 
other epochs in conjunction with bright thermal emission 
(a > 0.5) in the extension seen to the SE of the centrai 
source in the 8.5GHz map in Fig. ^[ This feature may be a 
result of locai condensation in the mass flow. 

In 1999, almost the entire emitting region has a nega- 
tive spectral index, with a minimum of a = —0.75 ± 0.07. 
This observation was made during what may be an eclipse of 
the inner symbiotic region by the outer giant (Skopal 1997). 
This may have reduced the relative intensity of the ther- 
mal emission, suggesting that the thermal emission region 
was very compact at this time, i.e. similar in extent to the 
angular diameter of the outer red giant. 



4 DISCUSSION 

4.1 The geometry of the extended emission 

The extensions detected in the radio surface brightness maps 
are more likely to be collimated jets rather than a torus. A 
torus is unlikely to appear so obviously one-sided, following 
an initial bipolar morphology. The non-thermal emission is 
suggestive of high velocity flow, commonly associated with 
jets but not with a circumbinary torus. 

It may be that the extended emission is due to a 
large nebula, with a radius of order 100AU, which has 
been photoionised between 1985 and 1986. From the STB 
model, for suitable parameters (Th = 100, 000K, M — 
1.9 x IO -6 Meyr" 1 (Skopal et al., 1996), K, ind =20 km s" 1 
(Vogel et al., 1994)) it is clear that the ionised cavity due 
to the hot component will be density bounded and much 
smaller than the binary separation. 



4.2 Expansion velocity 

The Hipparcos distance to CH Cygni is 268 ± 66 pc (Viotti 
et al. 1998), meaning that the maximum extent of the 3er rma 
contour of radio emission at 15GHz in 1986 is ~ 500AU. 
The most likely origin of of extended non-thermal emission 
is synchrotron emission from relativistic electrons in a mag- 
netic field. 

If CH Cygni follows the model proposed for HM Sagit- 
tae by Eyres et al. (1995), then what is being observed are 
regions of shocked material, caused when high velocity ejecta 
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Figure 3. VLA spectral index maps (centre) createti from 4.8GHZ (loft) and 8.5GHz or 15GHz (right) maps from from ali epochs. The 
lowest contour line on the brightness maps is cr TInB X —3, 3, 6, 12, 24, 48, 96 where present. The values of <r rms for each map are given in 
Tablc il. The white to black range on the spectral index map is —0.5 a ^ +1.5 and the contours represent changes in a of 0.1. 



© 2001 RAS, MNRAS 000,[lf|| 



The symbiotic star CH Cygni. I. Non-thermal bipolar jets 5 





1986 

i 

I 


■ 1 

ì ì ì 


1 I I ì 



I I I 


1988 

1 i 







ì 

ì 

ì 

i 


1995 

i 




I 





1999 


ì 

I I 
1 

1 


i 



-1.5 -1 -0.5 0.5 1 1.5 

Position along slice from SE to NW (arcseconds) 

Figure 4. Plots of instantaneous spectral index variations along 
the lines marked out in Fig. ^, from SE to NW, for each of the 
maps, relative to the bright peak in the higher frequency radio 
brightness map. 



expand into the relatively dense circumtrinary envelope cre- 
ated by the winds of the giant stars in the system. 

Evidence for this scenario comes from infrared spec- 
troscopy (Taranova & Yudin, 1988), carried out around the 
period of the end of the optical outburst in 1983-1987, show- 
ing changes in the circumtrinary dust envelope. They sug- 
gested that the dust was disrupted at the time of the radio 
jet emergence, probably as it was swept up by the expanding 
matter. 

The positions of the NW and SE extended knots were 
measured on the 1986 15GHz radio map. The angular sepa- 
rations between the peak in each knot and the centrai bright 
peak were 

A6W = 0.70 ± 0.01 arcsec 

and 

A0SE = 0.71 ± 0.05 arcsec 

for the NW and SE knots respectively. The greater uncer- 
tainty in the SE position is due to the fact that the SE com- 
ponent has lower surface brightness than the the NW one. 
In a 1985 VLA 15GHz map (Fig. these two knots were 
not visible, and hence the material must have been ejected 
within 422 days before the 1986 observation. This places a 
lower limit on the expansion velocity. Assuming a distance 
to CH Cyg of 268 ± 66 pc the velocities are found to be 

Vnw > 1210 ± 300 km s" 1 
for the NW component and 

Vse > 1230 ± 320 km s" 1 

for the SE component, with the uncertainties dominated by 
the uncertainty on the Hipparcos distance. These results 
correspond well with those derived previously, in particu- 
lar with the expansion velocity of 1400 km s _1 derived from 
two earlier maps (Taylor et al. 1986) that were separated by 
only 75 days. 

4.3 Magnetic fìeld in 1986 

The upper limit of the surface magnetic fìeld of a typical 
white dwarf is around 10G (Bond and Chanmugam, 1982). 
Assuming a dipole geometry, this will fall off as Bwd oc 
r~ 3 . For a stellar radius of O.OIRq the upper limit on the 
magnetic fìeld at 500 AU is 80 fiG, which is far too low to 
give efncient synchrotron emission. 

The speed of the unshocked pre-existing cool wind is 
likely to be much less than the jet speed, as a typical cool 
giant wind has a velocity of the order of 20-40 km s _1 (Vo- 
gel et al. 1994). The non-thermal emission is then caused by 
particle acceleration at the resulting bow shock, as the locai 
magnetic fìeld is compressed and enhanced and electrons are 
accelerated to velocities close to c. In such an environment, 
the locai azimuthal fìeld, H<f, is defìned by Equation 42 of 
Bode & Kahn (1985) the density of the shock, the tempera- 
ture of the giant wind (assumed to be Constant at ~ 10 4 K) 
and the mean particle mass which is taken to be the proton 
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Figure 5. VLA 15GHz map from 1985 January 22. The contours 
are as in Fig|Ì| and the rms noise level is 0.19mjy per beam. 



lie almost in the piane of the sky (Crocker et al., in prepa- 
ration), (f> ~ 90°. The cutoff frequencies of the radio emis- 
sion, z^mìn and ^ ma x are taken to be 0.01GHz and 100GHz 
respectively (Miley, 1980). The value of K is the main un- 
known in this equation as the exact details of the mechanism 
leading to synchrotron emission are not known. A value of 
K — 100 would be appropriate for the electrons being pro- 
duced following collisions in a circumstellar medium (Pa- 
cholczyk 1970), but a value of K = 1 is more consistent 
with the minimum energy requirement (Miley 1980). How- 

2 

ever, since H only depends upon K ? , having K 100 times 
larger would only increase He by a factor of 3. 

In order to determine the size of the jet bullets, two 
dimensionai gaussians were fitted to the emitting regions in 
the NW and SE jets. The depth was then found by assuming 
that the regions exhibited cylindrical symmetry along the 
axis of ejection. For the NW jet, 8\ — 0.39 arcsec, 82 = 0.34 
arcsec and d = 91 AU. The SE jet had dimensions of 81 — 
0.37 arcsec, 82 = 0.35 arcsec and d = 94 AU. Using the mean 
value of the spectral index in the jets of a — —0.23 ± 0.03, 



mass. The shock is assumed to be strong and totally ionised, 
so that the shock density is equal to 4 times the density of 
the unshocked wind at this distance from the star. It is as- 
sumed that the disruptions caused to the giant wind by the 
hot component prior to mass ejection are negligible at this 
distance so that a uniform, spherical radiai wind exists. 

The denisty of the wind is proportional to the mass loss 
rate of the giant, M = 1.9 x 10" 6 M Q yr" 1 (Skopal, 1996). 
If the unshocked wind speed, v w is taken to be the escape 
velocity of the red giant, assuming Af g i ant = O.93M0 and 
iìgiant = 123fl Q (Skopal 1997) then v w = 53.8 km s' 1 and 
H4 ~ 2.4 mG. 

Another magnetic field estimate can be derived from 
the technique used in studying high energy extragalactic 
sources. The total energy of a radio source is split between 
that of the electrons, the protons and the magnetic field. At 
equipartition, the magnetic energy is approximately equal 
to the particle energy (Pachoczyk, 1970), and the ratio of 
energy carried by protons to that carried by electrons is K. 
For a radio source of known angular dimensions and depth, 
assuming an ellipsoidal geometry, the magnetic field under 
equipartition, He, can be found using the method described 
by Miley (1980). 

To make use of this model, one needs to determine sev- 
eral parameters that describe the region: filling factor, 77; the 
angle between the magnetic field and the line of sight, \\ the 
flux density of the emitting region, S v (Jy), at frequency v 
(GHz) and the maximum and minimum frequencies of ra- 
dio emission (i/ max and Vmin repectively, in GHz) from the 
source spectrum. 

The filling factor would be more useful in extragalactic 
work when it would be expected that there would be unre- 
solved blobs of emitting material. Since the measured sizes 
of the knots in the extended emission are of a similar angu- 
lar size to the beam, 77 = 1. The angle \ is important due 
to the fact that the radiation observed is beamed from elec- 
trons moving towards the observer. The fact that eclipses are 
seen in the system means that the orbitai piane of the stars 
is observed edge-on. Bipolar jets are believed to be aligned 
with the rotation axes of their parent star (perpendicular to 
the orbitai piane) so the jets of CH Cygni are believed to 



H E = 3.8 mG 

for the NW jet and 

He = 3.6 mG 

in the SE. These values are comparable to that found above 
from the physics of the shock, so there is suffìcient magnetic 
field in the enhanced regions for synchrotron radiation. 

4.4 Lifetimes 

Since the total luminosity of this region at ali frequencies is 
not known, determining the lifetime of the emission due to 
ali radiative energy losses is impossible. In addition to this, 
synchrotron loss processes may be dominant. Under the in- 
fluence of a magnetic field H , an electron rotates at the gy- 
rofrequency (Rybicki & Lightman, 1979) and the frequency 
of radiation, v, is a factor of 7 3 higher than the gyrofre- 
quency, so manipulating Equation 6.4 of Rybicki & Light- 
man (1979) gives 7 = 1148. An electron with this Lorentz 
factor has energy r ym e c 2 and radiates with power P so that 
the lifetime of the synchrotron radiation is 

9m 3 e c 5 

and with the values of 7 and H = H$ from above, r = 
4.7 x IO 10 seconds, corresponding to just under 1500 years. 
Sucri a timescale would imply that synchrotron radiation 
would be very long-lived in the circumstellar material, as- 
suming that the magnetic field compression is sustained by 
continuous matter injection. When there is no matter in- 
jection the magnetic field can expand adiabatically and the 
synchrotron emission will fade within years. No non-thermal 
emission is seen in 1988. If indeed there is no synchrotron 
emission at this epoch, some other loss-mechanism must be 
occurring. The optical light curve of figure [j] shows 1988 to 
be a particularly quiet period, over two years after the last 
major outburst. It is therefore likely that there is no matter 
being ejected into the circumstellar wind, so that no new 
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density enhancements and magnetic field compressions are 
being formed. 

Supporting evidence for this scenario comes from opti- 
cal spectral observations. During this epodi, the blue con- 
tinuum was very faint and the Balmer lines were not present 
(Bode et al. 1991), strongly implying that there were also 
no fast winds or jets at this time. Hence one would not ex- 
pect to detect shock phenomena resulting from penetration 
of high-velocity material into circumbinary material and, by 
this model, no non-thermal emission from recently ejected 
material. 



4.5 Comparison with other symbiotic stars 

Previously, unambiguous non-thermal emission in symbiotic 
stars had only been seen in HM Sge (Richards et al. 1999). 
Spatially resolved non-thermal emission is very rare in any 
stellar source, and is only seen regularly in very high-energy 
systems such as micro-quasars. A potential non-thermal fea- 
ture in R Aqr was subject to such large errors (a a ~ 1) that 
its status could not be confirmed (Dougherty et al. 1995). 
Other symbiotic stars are believed to exhibit entirely ther- 
mal radiation. CH Cyg displays the strongest evidence yet of 
non-thermal emission due to jet-wind interaction. CH Cyg 
and HM Sge share the property of having undergone out- 
bursts in relatively recent history. The above arguments sug- 
gest that non-thermal emission will decrease rapidly with 
time after the outburst if the field strength and high energy 
particles are not continuously replenished as discussed ini- 
tially in Richards et al. (1999). Since the 1984-86 ejections 
in CH Cygni two outbursts have been observed which may 
have injected further high velocity material into the circum- 
stellar environment and hence given rise to the non-thermal 
emission seen in the 1995 and 1999 observations. 



5 CONCLUSION 

There is definite evidence of non-thermal emission in the 
ejecta from CH Cygni. This, together with a morphology 
consistent with high velocity jets, supports the scenario of 
Taylor et al. (1986), rather than one where the extension is 
a result of an expanding circumstellar torus. The proposed 
mechanism for the creation of the non-thermal emission in- 
volves material ejected in high velocity bullets interacting 
with a previously expelled lower velocity stellar wind. The 
ejected matter creates shocks in the circumstellar medium, 
leading to locally enhanced magnetic field density, particle 
acceleration and synchrotron emission. We suggest that the 
presence of non-thermal emission in symbiotic stars may be 
an indicator of outbursts associated with significant ejecta. 
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